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Abstract 

Chemical abundance of r-process elements in nearby dwarf spheroidal (dSph) galaxies is a powerful tool to 
probe the site of r-process since their small-mass scale can sort out individual events producing r-process elements. 
A merger of binary neutron stars is a promising candidate of this site. In faint, or less massive dSph galaxies such 
as the Draco, a few binary neutron star mergers are expected to have occurred at most over the whole past. We 
have measured chemical abundances including Eu and Ba of three red giants in the Draco dSph by Subaru/HDS 
observation. The Eu detection for one star with [Fe/H]=—1.45 confirms a broadly constant [Eu/H] of ~ —1.3 for 
stars with [Fe/H] xL —2. This feature is shared by other dSphs with similar masses, i.e., the Sculptor and the Carina, 
and suggests that neutron star merger is the origin of r-process elements in terms of its rarity. In addition, two very 
metal-poor stars with [Fe/H]=—2.12 and —2.51 are found to exhibit very low Eu abundances such as [Eu/H]< —2 
with an implication of a sudden increase of Eu abundance by more than 0.7 dex at [Fe/H] « —2.2 in the Draco dSph. 
In addition, the detection of Ba abundances for these stars suggests that the r-process enrichment initiated no later 
than the time when only a few % of stars in the present-day Draco dSph was formed. Though an identification of 
the origin of an early Eu production inside the Draco dSph should be awaited until more abundance data of stars 
with [Fe/H] ;$ —2 in the Draco as well as other faint dSphs become available, the implied early emergence of Eu 
production event might be reconciled with the presence of extremely metal-poor stars enriched by r-process elements 
in the Galactic halo. 
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1. Introduction 

Neutron star (NS) mergers as the origin of r-process ele¬ 
ments have been highlighted thanks to the recent updates from 
both observational and theoretical studies. First of all, the de¬ 
tection of a bump of a near infrared light in the afterglow of a 
short-duration 7 -ray burst has become compelling evidence for 
the r-process synthesis in a NS merger event (Barnes & Kasen 
2013; Tanvir et al. 2013; Berger et al. 2013; Hotokezaka et al. 
2013). Another advancement of our knowledge is brought by 
the results of theoretical nucleosynthesis calculations; an alter¬ 
native candidate for the r-process site, neutrino-driven winds 
of core-collapse supernovae (CCSNe), have been found to be¬ 
come proton-rich due to revised treatments of neutrino transfer 
and unable to synthesize heavier r-process elements (A>130) 
(Thompson et al. 2001; Wanajo 2013) while recent studies have 
shown that NS mergers synthesize and eject them (Korobkin 
et al. 2012; Bauswein et al. 2013; Wanajo et al. 2014). In 
addition, chemical feature of r-process elements (Eu and Ba) 
in the Galaxy can be reproduced by enrichment proceeded 
through NS mergers in a novel model (Tsujimoto & Shigeyama 
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2014, hereafter TT14), which improves the mixing process of 
r-process elements from the previous ones (Mathews & Cowan 
1990; Argast et al. 2004) to reflect the kinematic features of 
ejecta of NS mergers. 

TT14 have proposed a more direct assessment of r-process 
origin by distinguishing between CCSNe and NS mergers from 
observed chemical abundances in nearby dwarf spheroidal 
(dSph) galaxies. Since NS merger events are considered to be 
rare, galaxies with small masses such as faint dSph galaxies are 
unlikely to leave so many NSs after the deaths of massive stars 
as to realize their mergers. For instance, a galaxy with its stel¬ 
lar mass of 10 5 M 0 are expected to host ~ 500 CCSNe in total 
assuming a canonical stellar initial mass function (IMF). Then, 
the estimated frequency of NS mergers of one per 1000-2000 
CCSNe (TT14) <1 NS merger in this galaxy. TT14 exam¬ 
ined the observed trend of [Eu/H] against [Fe/H] for three faint 
dSphs, that is, the Draco, the Carina, the Sculptor, and found 
a constant [Eu/H] of ~ —1.3 with no apparent increase in the 
Eu abundance over the metallicity range of — 2 [Fe/H] /$ — 1. 
This result suggests that the r-process production is associated 
with events much rarer than SNe, such as NS mergers. On the 
other hand, bright (i.e., massive) dSphs such as the Fornax and 
the Sagittarius exhibit a clear increasing [Eu/H] trend with re¬ 
spect to [Fe/H] as in the Milky Way, as might be expected since 
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Name 

RA 

DEC 

V 

Date 

ExpTime (sec) 

S/N 1 

Vi os (km s L ) 

Draco targets 

Irwin 19826 

17 20 52.97 

+57 55 57.9 

17.23 

2014/08/16 

12900 

52 

-292.5 + 0.7 

Irwin 20751 

17 22 13.63 

+57 53 06.5 

17.39 

2014/08/18 

16200 

48 

-301.5 + 0.5 

Irwin 21275 

17 19 41.84 

+57 52 19.2 

16.87 

2014/08/17 

15300 

66 

-295.0 + 0.4 

Comparison halo stars 

BD+30 2611 

15 06 53.83 

+30 00 36.94 

9.13 

2014/08/18 

600 

238 

-281.7 + 0.3 

HD 122956 

14 05 13.02 

-14 51 25.45 

7.25 

2014/08/16 

200 

195 

165.3 + 0.3 

BD+09 2870 

14 16 29.95 

+08 27 52.94 

9.70 

2014/08/17 

600 

325 

-119.7 + 0.2 


i a signal-to-noise per pixel (~ 5 pixels per resolution element) measured by taking the standard deviation in 
the wavelength range 6649-6653A. 


Table 1 . Summary of the observation 


NS mergers are not so rare in these bright galaxies. 

One puzzling question on this dSph matter is why and how 
these dSph stars are so promptly enriched with r-process el¬ 
ements up to [Eu/H]~ —1.3, while these stars claim that no 
r-process productions operate during a major star formation 
epoch with [Fe/H] >, —2 (see $ 4). The idea that the gas form¬ 
ing these galaxies was already enriched at a level of [Eu/H |~ 
— 1.3 should be discarded since very metal-poor stars exhibit 
lower Eu abundance as observed for the star Draco 119 show¬ 
ing [Eu/H]< -2.55 with [Fe/H]=-2.95 (Fulbrightet al. 2004). 
A few stars with [Fe/H]< —3.5 in the Sculptor also exhibit 
[Eu/H]< -2.2 (Tafelmeyer et al. 2010; Simon et al. 2015). 
Such low abundances of r-process element is verified by the 
Ba abundance for very metal-poor ([Fe/H]< —2) dSph stars 
(Roederer 2013). Therefore it is likely that some kind of r- 
process synthesis may operate at a very early stage of star for¬ 
mation ([Fe/H] — 2) inside these dSphs and lift Eu abun¬ 

dance up to [Eu/H |~ —1.3. Irrespective of its production site, 
we need to answer how events of r-process synthesis that oc¬ 
curred at a very early stage could leave no trace in stars with 
[Fe/H] >, —2 which constitutes a main body of dSphs. 

We aim at making a thorough investigation of the trend of r- 
process elements imprinted in stars covering the whole metal- 
licity range in faint dSphs by measuring stellar Eu abundance. 
Note that Eu is an ideal element to trace the r-process enrich¬ 
ment; 97% of the solar Eu abundance was synthesized by r- 
process, while the fraction of r-process is only 15% in the solar 
Ba abundance (Burris et al. 2000). In particular, we will fo¬ 
cus on (i) confirming a constant [Eu/H] ratio of ~ —1.3 for 
[Fe/H] ^ — 2 and (ii) understanding how the [Eu/H] ratio in¬ 
creases up to ~ —1.3 for [Fe/H] + —2, with a sufficient number 
of stars in the Draco and the Sextans dSphs. 

This first paper presents results of chemical abundances 
measurements for three red giants in the Draco dSph and dis¬ 
cusses its theoretical implications for the origin and enrichment 
history of r-process elements. 

2. Observation 

Three target red giant stars belonging to the Draco dSph 
were selected from Kirby et al. (2011). For one of our target 
stars, Irwin 21275, an upper limit of Eu abundance has been 
estimated to be [Eu/H]< —1.8 (Cohen & Huang 2009), while 
no Eu measurements were available for the other two stars. 



6644 6646 6648 

Restframe wavelength A 


Fig. 1 . Spectra of the Draco stars and one of the comparison stars 
(BD+30 2611) including the Eu II 6645 A absorption line. The red 
lines show synthetic spectra that give the best fits to the Eu II line for 
Irwin 19826 (top panel) and BD+30 2611 (bottom panel). For other 
two stars (two middle panels), synthetic spectra con'esponding to a 20 - 
upper limit of the Eu abundance are shown. 

Observations were made with High-Dispersion Spectrograph 
(HDS) mounted on the Subaru Telescope (Noguchi et al. 2002). 
We adopted one of the standard setups of HDS, gStdYd”, 
which covers a wavelength range of ~ 4000 — 6800 A. We 
use a slit width of l."l to achieve a spectral resolving power 
of R~ 30000. The 2 x 2 binning was applied for our spectra. 
Summary of our observations of the three Draco stars together 
with comparison halo stars is tabulated in Table 1. 

Data were reduced with the HDS reduction pipeline ghd- 
sql”, which implements IRAF 2 scripts for cosmic-ray removal, 

2 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National Science 
Foundation. 
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No.] 


Name 

T e fi 

log g 

5 

[Fe/H ] 3 

log eEu 

[Eu/H ] 4 

log eBa 

[Ba/H] s 


(K) 

(dex) 

(km s 4 ) 

(dex) 

(dex) 

(dex) 

(dex) 

(dex) 

Irwin 19826 

4028 

0.6 

1.7 

— 1.45 ± 0.12 

-0.81 

-1.33 ±0.21 

0.04 

—2.14 ± 0.23 

Irwin 20751 

4278 

0.6 

1.7 

—2.12 ± 0.13 

< -1.5 

< - 2.0 

-0.41 

-2.59 ±0.15 

Irwin 21275 

4406 

0.7 

2.3 

-2.51 ±0.09 

< -1.6 

< - 2.1 

-0.99 

—3.17 ± 0.13 

BD+30 2611 

4201 

0.6 

1.8 

— 1.48 ± 0.12 

-0.47 

-0.99 ±0.14 

- 

- 

HD 122956 

4595 

1.3 

1.7 

— 1.72 ± 0.12 

- 0.88 

-1.40 ±0.14 

0.52 

— 1.66 ± 0.16 

BD+09 2870 

4612 

1.2 

2.0 

-2.47 ±0.09 

< —2.45(—2.64) 

< —2.97(—3.16) 

-0.93 

—3.11 ± 0.11 


3 The weighted average of the estimates from neutral ([Fe I/H]) and ionized ([Fe II/H]) species. The uncertainty in each of 
the [Fe I/H] and [Fe II/H] values is calculated as a quadratic sum of contributions from line-to-line scatter (root-mean 
square difference from the mean value divided by a square-root of the number of lines) and uncertainties in atmospheric 
parameters, which are assumed to be A T e s = 100 K, Alogg = 0.3 dex, and A£ = 0.2 km s -1 . 

4 The values in the parenthesis are deduced from the Eu II 4129.7 A line, while other values are from Eu II 6645.1 A. 
Adopted solar abundance is logeEuQ = 0.52 (Asplund et al. 2009). The uncertainty includes contributions from a 
statistical error and uncertainties in the atmospheric parameters. 

5 Adopted solar abundance is logeBaQ = 2.18 (Asplund et al. 2009). The uncertainty is calculated by the same way as for 
[Fe I/H]. 


Table 2. Atmospheric parameters and derived abundances 


flat fielding, scattered-light subtraction, and wavelength cali¬ 
bration. The reduced spectra are normalized by fitting the con¬ 
tinuum fluxes with polynomial functions. Examples of the re¬ 
duced spectra of the three program stars together with one com¬ 
parison star are shown in Figure 1. 

3. Abundance analysis 

We performed a standard abundance analysis using a 1D- 
LTE analysis code adopted in Aoki et al. (2009), which imple¬ 
ments model atmospheres of Castelli & Kurucz (2003). The 
procedure is briefly summarized as follows (see Ishigaki et al. 
2014 for details). First we measure equivalent widths of ab¬ 
sorption lines by fitting a Gaussian profile to each absorption 
feature. The effective temperatures are estimated by (V — K) 
color and an assumed Fe abundance, using the empirical for¬ 
mula for red giants (Ramirez & Melendez 2005). The log g 
values for the Draco stars are determined using the estimated 
T e s by assuming that each star is located on a red giant branch 
of Yonsei-Yale isochrones (Demarque et al. 2004) with an age 
10.9 Gyr (Orban et al. 2008). Then we deduce Fe abundances 
from the equivalent widths of individual Fe I lines. Here the 
micro-turbulence velocity parameter (0 is determined so as to 
meet the condition that Fe abundances dedudced from Fe I lines 
with different excitation energies do not depend on the equiva¬ 
lent widths. With these stellar parameters, we revise the effec¬ 
tive temperature and iterate the same procedure until a consis¬ 
tent set of stellar parameters including the Fe abundance is ob¬ 
tained. The results are summarized in Table 2. It should be of 
note that the estimated parameters for Irwin 21275 agree well 
with those estimated by Cohen & Huang (2009). The differ¬ 
ences between the two analyses are AT e g = 90 K, Alogg = 0.0 
dex, A£ = 0.1 km s _1 , and A [Fe/H]= 0.1 dex. 

The Eu abundance is estimated by fitting synthetic spectra 
to each absorption feature. The wavelength region around the 
Eu II 4129.72 A line is too noisy for the Draco stars to define 
a continuum level. Instead we use the Eu II 6645.10 line for 
the Eu abundance estimate. The local continuum level is de¬ 


fined from the wavelength range of 6649.0-6653.0 A, in which 
strong absorption lines are expected to be absent for the stellar 
parameter range of our target stars. We inspected CN features 
in the spectra of the three stars, and found that their contribu¬ 
tion is negligible in this wavelength range. The modification 
of the line profile due to the hyperfine structure is taken into 
account assuming an isotopic fraction for the two naturally ex¬ 
isting isotopes, 151 Eu and 153 Eu, as 151 Eu/( 151 Eu± 153 Eu) = 
0.5 (Lawler et al. 2001). Results of the fitting are shown in 
Figure 1. As shown in the two middle panels, the Eu II lines 
are too weak to be detected for the two stars with [Fe/H] < —2.0 
and thus only upper limits are deduced. On the other hand, for 
one metal-rich star with [Fe/H] = —1.45, the equivalent width 
due to the Eu II feature is estimated to be 20 mA while the 
2 cr uncertainty estimated from the noise level of the contin¬ 
uum wavelength range is approximately 6 mA. Thus we ar¬ 
gue that we detected the Eu II line and that its abundance is 
[Eu/H]= —1.33. The uncertainty in the equivalent width corre¬ 
sponds to an uncertainty in the Eu abundance of no more than 
0.15 dex. Changing each of the adopted atmospheric param¬ 
eters by AT e ff ± 100 K, Alogg ± 0.3 dex, and A£ ± 0.2 km 
s _1 results in A[Eu/H] of 0.13 dex at most. The Ba abun¬ 
dances are estimated from the measured equivalent widths of 
Ba II 5853.69, 6141.73, 6496.91 A absorption lines. The ef¬ 
fect of the hyperfine structure is also taken into account in the 
same way as in Ishigaki et al. (2013). Results of the abundance 
measurements are summarized in Table 2. Detailed elemental 
abundances including a-, Fe-group, and other neutron-capture 
elements will be presented in our forthcoming paper (Ishigaki 
et al., in preparation). 

4. Eu feature and its implication 

In this section, we present the chemical feature of Eu abun¬ 
dance revealed by this study together with previous observa¬ 
tions for not only the Draco dSph but also other two dSphs, 
i.e., the Sculptor and the Carina. Then, we try to advance our 
knowledge on the enrichment history of r-process elements as 
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Fig. 2. Observed Eu abundances against Fe abundances in the Draco 
dSph. Our results are indicted by large filled circles, together with 
previous results by Shetrone et al. (2001), Fulbright et al. (2004), and 
Cohen & Huang (2009). 

well as their origin. First we show our results for the Draco 
dSph in Figure 2. The Eu feature can be summarized as (i) the 
broadly constant Eu abundance of [Eu/H]~ —1.3 for [Fe/H] 
^ — 2 and (ii) low Eu abundances of [Eu/H] /$ — 2 deviating 
from a constant value for [Fe/H] ;$ —2.3-2.1. A major¬ 

ity of stars in the Draco dSph show the first feature, i.e, the 
constant Eu abundance, because the stellar metallicity distribu¬ 
tion function (Kirby et al. 2011) indicates that ~ 85% of stars 
populate the metallicity range of [Fe/H]> —2.3 or ~ 70% of 
stars within [Fe/H]> —2.1. Thus the constant Eu abundance for 
[Fe/H] £ —2 suggests a rarity of r-process production events. 
Using the Kroupa IMF (Kroupa et al. 1993), we obtain the cor¬ 
responding stellar mass of ~ 2.5 x 1O 5 M 0 (Martin et al. 2008) 
and find that this stellar population presumably hosted ~1400 
CCSNe in total. Accordingly we conclude that the r-process 
production did not happen during sequential events of ~ 1400 
CCSNe in the Draco dSph. This implied rate is in good agree¬ 
ment with that of NS mergers, one per 1000 — 2000 CCSNe, 
which is deduced from an analysis of chemical abundances in 
massive dSph galaxies (the Fornax and the Sagittarius) as well 
as in the Milky Way (TT14). 

This broadly constant Eu feature breaks at a very low- 
metallicity range. Our result for Irwin 20751 indicates a very 
low [Eu/H] by more than ~ 0.7 dex at [Fe/H]=—2.1 ± 0.13 as 
compared with the plateau abundance of [Eu/H]~ —1.3. Since 
the star with [Fe/H]=-2.34± 0.19 (Cohen & Huang 2009) al¬ 
ready belongs to the plateau-like Eu abundance branch, [Eu/H] 
is likely to have an abrupt change at [Fe/H]« —2.2. Other 
two stars with [Fe/H]< —2.5 exhibiting low Eu abundances 
clearly confirm the presence of another branch of [Eu/H] for 
[Fe/H] ;$ —2 far deviating from a constant [Eu/H] of —1.3. 
Note that Ba is detected in the spectra of our two stars with 
[Fe/H]< —2 as well as the star with [Fe/H]=—2.95 (Fulbright 
et al. 2004). This implies that these stars contain some amount 
of Eu. Accordingly, from an overall feature of Eu abun¬ 
dance, we conclude that the r-process operates inside the Draco 
dSph only at an early star formation with its enrichment up to 
[Eu/H]~ —1.3, while a majority of stars formed in the subse¬ 
quent epoch did not yield r-process elements. Regarding the 



Fig. 3. observed Eu feature for the Draco (red circles). Sculptor (blue 
circles), and Carina (crosses) dSphs. The results acquired by this study 
are indicted by big red circles. Observed data are taken from Shetrone 
et al. (2003); Geisler et al. (2005); Kirby & Cohen (2012); Simon et al. 
(2015) for the Sculptor dSph, and Shetrone et al. (2003); Lemasle et al. 
(2012); Venn et al. (2012) for the Carina dSph. 

site of an early Eu production, the feature of a sudden increase 
in [Eu/H] by more than 0.7 dex associated with little increase in 
[Fe/H] occurring around [Fe/H]=—2.2 seems compatible with 
the enrichment by a NS merger. 

A common feature is shared with the Sculptor and Carina 
dSphs (Fig. 3). Clearly we see a constant [Eu/H] for 
[Fe/H] —2 for the three dSphs. For [Fe/H] ;$ —2, the Sclptor 
dSph also exhibits the low Eu abundance deviating from a 
constant value by more than 3 dex, though the feature of a 
gradual increase in [Eu/H] between [Fe/H]« —2.5 and —2 
is different from that of the Draco dSph. In addition, we 
present the evolution of [Ba/Eu] for the three dSphs in Figure 
4. We see that dSph stars exhibit a clear increasing trend of 
[Ba/Eu], which is considered to be a result of a gradual en¬ 
richment of Ba by s-process, though a few stars including our 
star with [Fe/H]=—1.45 deviate from this trend. Since the ob¬ 
served lower limits of [Ba/Eu] ratios for the two stars with 
[Fe/H]< —3 are consistent with a pure r-process Ba/Eu ratio of 
[Ba/Eu]=—0.89 (Burris et al. 2000), we may be able to guess 
the values of [Eu/H] from their Ba abundances. Assuming that 
these two stars with [Fe/H]=—4.05 and —3.52 hold the pure r- 
process Ba/Eu ratio, we deduce individual [Eu/H] abundances 
of —4.57 and —3.14, respectively. It may give the bottom level 
of [Eu/H] abundance in faint dSphs. 

Still too small data for [Fe/H] ;$ —2.5 cast a veil over the 
evolutionary path of Eu abundance from the bottom, likely 
[Eu/H]< —3, to [Eu/H]~ —1.3 and prevent us from identifying 
the site of the early Eu production inside faint dSphs. Even if 
we add information on the Ba feature (Fig. 5), its clear answer 
continues to elude us. It seems questionable to assume a selec¬ 
tive early NS merger which should occur in each of the three 
dSphs, given its rarity. We may need to introduce some other 
contributors as an early r-process production such as magneto- 
rotationally driven SNe (Winteler et al. 2012; Wehmeyer et al. 
2015). To advance our understanding of the origin of early r- 
process production, more Eu data with precise measurements 
for stars with [Fe/H] — 2 in faint dSphs together with theo- 
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Fig. 4. Observed [Ba/Eu] ratios for the three dSphs. Symbols are the 

same as those in Fig. 3. 

retical updates of the models for both NS merger and SNe are 
surely awaited. In particular, we can raise future prospect that 
improving the S/N values such as >, 30 at the blue wave length 
of Eu II 4129.7 A likely achieved by 30m-class telescopes will 
reveal the Eu feature for very metal-poor regime. Putting aside 
an enigma on the site of an early Eu production, a clear signa¬ 
ture of Eu enrichment imprinted on very early generations of 
dSph stars is reconciled with the fact that extremely-metal-poor 
halos stars in the Galaxy are already enriched in r-process, im¬ 
plying little time delay of r-process ejection as compared with 
other elements such as a-elements or Fe. 

5. Concluding remarks 

From an assembly of Eu abundance data thus far obtained 
for three faint dSph galaxies including our new data for the 
Draco, it is clear that the trend of [Eu/H] against [Fe/H] shows 
completely different behaviors between two episodes, i.e., a 
very early phase and the subsequent major star formation pe¬ 
riod; a remarkable evolution of [Eu/H] from <—4to~—1.3 
in the metallicity range of [Fe/H] A 2 is followed by broadly 
constant [Eu/H] until the end of star formation. The latter fea¬ 
ture exhibited by a main body of dSph stars that host more 
than 1000 CCSN occurrences implies an extremely low event 
of the Eu production. Though NS merger is thought to satisfy 
this rarity as a source of r-process elements, there exist the Eu 
production events in the early phase of dSph formation, which 
remains a mystery. Such an early commencement of r-process 
enrichment is likely to have something in common with the 
early r-process enrichment of the Galactic halo. 

We are grateful to A. Tajitsu, K. Aoki, and staff members 
of the Subaru telescope for their helpful support and assistance 
in our HDS observation. This work was supported by JSPS 
KAKENHI Grant Number 23224004, 25-7047. 
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